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The conductance, the transmission and the reflection probabilities through rectangular potential 
barriers and pn-junctions are obtained for bilayer graphene taking into account the four bands of 
the energy spectrum. We have evaluated the importance of the skew hopping parameters 73 and 
74 to these properties and show that for energies E > 71/IOO their effect is negligible. For high 
energies two modes of propagation exist and we investigate scattering between these modes. For 
perpendicular incidence both propagation modes are decoupled and scattering between them is 
forbidden. This extends the concept of pseudospin as defined within the two band approximation 
to a four band model and corresponds to the (anti)symmetry of the wavefunctions under in-plane 
mirroring. New transmission resonances are found that appear as sharp peaks in the conductance 
which are absent in the two band approximation. The application of an interlayer bias to the 
system: 1) breaks the pseudospin structure, 2) opens a bandgap that results in a distinct feature of 
suppressed transmission in the conductance, and 3) breaks the angular symmetry with respect to 
normal incidence in the transmission and reflection. 

PACS numbers: 73.22.Pr, 72.80.Vp, 73.63.-b 



I. INTRODUCTION 

Bilayer graphene (BLG) is a system consisting of two 
BernaP stacked graphene monolayers^. Whereas mono- 
layer graphene has a linear electronic spectrum in the 
vicinity of each corner of the Brillouin zone (K and K' 
points, also known as Dirac points), bilayer graphene has 
four hyperbolic bands. Two of these bands touch in the 
K point at zero energy making BLG a gapless semicon- 
ductor, the other two bands are displaced by an energy 
of 71 = 377meV with respect to the touching bands . 
The application of a potential that breaks the interlayer 
symmetry can however create a tunable bandgapP^. Of- 
ten, a low energy approximation^ is made that is valid for 
electron kinetic energy much smaller than the interlayer 
hopping parameter 71 . This so called two band approxi- 
mation has a quadratic dispersion and is only valid near 
the Dirac point. 

In monolayer graphene Klein tunneling results in 
a 100% probability for perpendicular transmission 
through pote ntial b arriers, as predicted^ and observed 
experimentally^^. For bilayer graphene due to the con- 
servation of pseudospin, no Klein tunneling is expected 
and this was confirmed theoretically within the two band 
approximation^!^. i n this case there are electronic states 
available that are not accessible for penetration into the 
potential barrier, which was called a cloaking of those 
states^. 

Previous work was based on the two band approxi- 
mation which we extend here to the four band model. 
This allows us to investigate the electronic properties at 
higher Fermi level, i.e. beyond 71, and for a higher elec- 
trostatic potential. Recent experimental progress has al- 
lowed to access this energy region and measurements of 
the electronic transport in this region is expecte d 14 " 16 !. 
We calculate the transmission and reflection probabil- 
ities within the same band and between the two bands 



for electrons impinging on a rectangular potential barrier 
(pnp-junction) and a potential step (pn-junction) at dif- 
ferent angles of incidence and investigate the effect of the 
application of an interlayer bias to the transport proper- 
ties. Furthermore we compare the energy dependence of 
the conductance calculated within the two band and the 
four band model and point out major differences. 

Our findings show that at low energy and low poten- 
tial the same phenomena occur as those predicted by the 
two band model and that the existence of pseudospin is 
related to the (anti) symmetry of the wavefunctions. This 
relation is a consequence of the symmetry of the crystal 
and is therefore also valid when the skew hopping pa- 
rameters arc taken into account. Outside the range of 
validity of the two band approximation however, we pre- 
dict several new phenomena: 1) for high potential bar- 
riers new resonances are found that are absent in the 
two band model; 2) the use of four bands introduces a 
new mode of propagation; 3) we investigate a new form 
of cloaking and calculate the scattering between the two 
modes of propagation; and 4) the newly discovered reso- 
nances and the second mode of propagation lead to dis- 
tinctive features in the conductance that are absent in 
the two band calculation. We also justify the use of only 
the nearest neighbor interlayer hopping parameters. Fi- 
nally, we show that the application of an interlayer bias 
not only opens a bandgap and therefore suppresses the 
transmission in this region, it also unexpectedly breaks 
the angular symmetry with respect to normal incidence 
when only one valley is considered. 

The paper at hand is organized as follows. In Sec. [TT] 
we present the formalism, indicate the different propa- 
gating modes, define the eight different transmission and 
reflection probabilities for the four band model and ex- 
plain the transition from the four band to the two band 
model at low energy. The effect of the skew hopping pa- 
rameters on our results is critically examined. In Sec. |III| 
we analyze the symmetries of the system to explain the 
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transmission for normal incidence and the surprising oc- 
currence of angular asymmetry at non normal incidence 
The numerical results for the conductance, transmission 
and reflection at non normal incidence for pn-junctions 
and potential barriers with and without interlayer bias 
are discussed in Sec. |IV[ and we show the effect of the 
skew hopping parameters on the transmission. In Sec. [V] 
we summarize the main points of this paper. 



II. THE PROPAGATING MODES 

In this section we discuss the dispersion relation of 
BLG and the resulting propagating modes. It turns out 
that electrons in BLG can propagate via two different 
modes and we find that when electrons impinge perpen- 
dicular on a potential barrier, it is not possible to scatter 
between those modes. 

We model the BLG crystal as two hexagonal monolayer 
flakes with in plane interatomic distance^ a = 0.142nm, 
each consisting of two nonequivalent sublattices with 
atoms Ai and B\ for the top layer and A2 and B2 for 
the bottom one. These two layers are stacked according 
to Bernal stacking which places the A 2 atom just above 
the B\ as schematically shown in Fig. [lja). In both lay- 
ers, each A atom is surrounded by three B atoms and 
vice versa. The intralayer coupling between these atoms 
is 70 ~ 3eV. Between the A 2 and Bi atoms the in- 
terlayer coupling is 71 ~ OAeV while the skew hopping 
energy between the other two sublattices are denoted as 
73 « 0.3eV and 74 ~ OAeV. These interatomic cou- 
pling parameters are depicted in Fig. [lja). The contri- 
bution of th e skew hopping parameter 73 results in the 
so callecpEsl trigonal warping, an effect occurring only at 
very low energy (E < 4meV). The other parameter, 74 
has an even lower impact on the electronic properties as 
will be discussed. Therefore, we will often neglect these 
two 7-parameters allowing for a more comprehensive dis- 
cussion. 

Following the continuum nearest-neighbor tight- 
binding formalism, the effective Hamiltonian near t he K 
point and the corresponding eigenstates are given b} * 19 * 20 ! 
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Here, vf = ~ 10 6 m/s is the Fermi velocity for elec- 
trons in monolayer graphene and 1)3,4 = 1^73,4/70 are 
related to the skew hopping parameters, ir = p x + ip y = 
fr(k x + iky) is the in-plane momentum relative to the 
Dirac point, V is a general potential term and 5 cor- 
responds to an externally induced interlayer potential 
difference. Due to the dimensionality of the Hamilto- 
nian, the eigenstate of the system is a four component 
spinor. Neglecting the skew hopping parameters, the en- 



ergy spectrum of this Hamiltonian is given by 
e = / I \ k 2 + A 2 
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We have used 
A = S/hvF and 



the reduced variables Ti = ^l/hvp 
e = (E — V)/Hvf with E the energy of the electrons. 
For a system without interlayer potential bias S, this re- 
sult reduces to the one previously found by Snyman et 
aZPP. The energy spectra corresponding to these systems 
are displayed in Fig. |ljc). In Fig. [ljb) we show the dis- 
persion relation when taking also account of the skew 
hopping parameters 73 and 74. Their effect is clearly 
negligible for E > 71/ 100 « 4meV. 

Vice versa it is possible to calculate the value of k x 
as a function of the energy and k y . This corresponds 
to the wave vectors of the plane wave solutions of the 
Schrodinger equation = E^ and is given by 



k l = We 2 + A 2 + iJe 2 (r 2 + 4A 2 ) - TfA 2 - k 2 (3) 



which for A = reduces to 



k l = Je 2 + leV l -k 2 v . 



(4) 



Depending on the value of £ relative to T\, A and k y this 
wave vector can be real or imaginary. This means that 
it can represent a traveling or an evanescent plane wave. 
Without the interlayer bias, when < £ < T%, the k~ is 
imaginary while k + is real and then propagation is only 
possible using the k + channel. When e > T\ however, k~ 
becomes real too, providing a new mode for propagation. 
For e < 0, a similar argument leads to propagation via 
k~ when |e| < and two ways of propagation when 

|£|>IY 

One obtains the two band approximation for e <C T\ 
with the condition that A and e are of the same order of 
magnitude. Therefore, one can neglect these terms up to 
second order in Eq. Q, resulting in the wave vector 
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which for A = reduces to 



k l « JleTi - k 2 , 
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FIG. 1. (Colour online) (a) Schematic presentation of the 
sublattices in bilayer graphene. The arrows indicate the dif- 
ferent interatomic hopping parameters, (b) Energy spectrum 
of bilayer graphene near one of the Dirac points for low en- 
ergy. The dotted curve corresponds to the spectrum account- 
ing for all interatomic hopping parameters, the solid curve 
accounts only for 70 and 71. (c) Total energy spectrum of 
bilayer graphene for (left) an unbiased system and (right) a 
system with interlayer bias 8 — O.371. The dotted curves 
also account for the skew hopping parameters while the solid 
curve considers only nearest neighbor interlayer hopping. The 
dashed curve corresponds to the spectrum of the two band ap- 
proximation. All bands in (b) and (c) are coloured according 
to their relation to the (anti)symmetric states. 



This spectrum is superimposed in Fig. |T](c) as dashed 
curves. It agrees with the full spectrum only near the 
Dirac point. The validity of the approximation is based 
on the increase in energy near the atomic sites of the Bi 
and Ai atoms which influence each other. For low Fermi 
energy, it therefore makes sense to take into account only 
the orbital wave functions near the other two atoms. This 



represent a traveling state when s < which corresponds 
to a hole state. In contrast to the four band treatment, 
there is no second mode of propagation. The introduction 
of an interlayer bias term S in the system has a strong 
influence on the electronic properties. As was pointed 
out earliei^, it opens a gap in the spectrum which com- 
pletely changes the low energy behavior of the electrons. 
When an interlayer bias is applied, the effect of the skew 
hopping parameters is of even less importance for the 
spectrum and as shown in Fig. [ljc). Notice that the two 
band approximation fails to describe the spectrum of the 
system accurately. 

The potential we will consider is similar to the one 
used by Katsnelson et alW, but now with the addition 
of an interlayer potential bias term. It consists of a one 
dimensional potential barrier of width d given by 

!0 if x < (region I) 

V +£6 if0<x<d (region II) (12) 
if x > d (region III) 

where £ = +1 for the first layer and £ = —1 for the 
second layer. This barrier is shown in Fig. [2] and be- 
cause it is translational invariant in the y direction, k y is 
a conserved quantity. Using the Hamiltonians given by 
Eqs. ([I]) and ([9]), we can now calculate the transmission 
and reflection probabilities for electrons impinging on the 
barrier. The two band model only allows for one mode 
of propagation, leading to one transmission (T) and one 
reflection (R) channel. For sufficiently large energy how- 
ever, in the four band model, it is possible to propagate 
through two distinct modes. Therefore, we will have four 
reflection and four transmission channels. For the trans- 
mission these are two non scattered channels, which we 
denote as and T~ for propagation via k + and fc~ re- 
spectively and two scattered channels in which the parti- 
cle enters through one channel and exits through another 
one. We will denote them as T_ for scattering from the 
k + band to the fc~band and for the other direction. 
A similar definition holds for the R+ reflection channels. 



reduces thejf x 4 Hamiltonian, Eq^ 0, and replaces it The eight channels are schematically depicted in Fig. [| 



with the approximate one given b 
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where V" = -71V/ (H 2 v 2 F ) and A" = - 7l <$/ (h 2 v 2 F ). 
The two spinor plane wave solution of the Schrodinger 
equation of this Hamiltonian consists of a propagating 
wave with wave vector given by Eq. Q and an evanes- 
cent mode with inverse decay length 



Using the transmission probabilities, we can calculate 
the conductance as function of the energy given by the 
Landauer-Biittiker formula^ 
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which corresponds to the imaginary part of k l with 
I = -1 in Eq. There is no positive energy value 

that can make this quantity imaginary and so it can only 



where Go = 4e 2 //i, which is four times the quantum of 
conductance due to spin and valley degeneracy and W is 
the width of the sample in the y direction. 



III. TRANSMISSION PROBABILITIES AND 
SYMMETRIES 

Before calculating the transmission probability by 
matching the components of the spinor wave function 
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tions Vl/^i which are described by the 2x2 Hamiltonian 
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FIG. 2. (Colour online) Schematic representation of the trans- 
mission and reflection probabilities for a rectangular potential 
barrier. 



at the boundary as explained by Barbier et aZP^ and 
Snyman et a/P^I, a simplification can be made by trans- 
forming the Hamiltonian from Eq. ([I]). Constructing 
symmetric and antisymmetric combinations of the spinor 
components by combining the atomic wave functions 
with ipB 2 and "0-Bi with ipA 2 , the Hamiltonian is trans- 
formed into 
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with the new spinor 
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where we have introduced the reduced potential V = 
V/Hvp and u 4 = (1 ± v' 4 ) with the dimensionless skew 
hopping velocities v 3 4 = V3. 4 / v f — 73,4/70 which turn 
out to be very small, i.e. v' 3 « 0.09 and v' 4 « 0.04. These 
small values advocate the neglect of the skew hopping pa- 
rameters in the discussion a little bit further. The four 
component spinor can be seen as a combination of two 
two-spinors ^4.+ and ^4,— which are respectively anti- 
symmetric and symmetric with respect to the exchange 
A\ f-> B 2 and A 2 f-> B 4 . This exchange corresponds to 
a reflection of the system by an in-plane mirror. Notice 
that for normal incidence and for an unbiased system, i. 
c. when k y = and A = 0, the Hamiltonian is block 
diagonal in this basis. This means that it represents two 
non interacting one dimensional systems with eigenfunc- 
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Neglecting the skew hopping terms, this Hamiltonian de- 
scribes a one dimensional monolayer of graphene with a 
potential term that breaks the sublattice symmetry and 
corresponds to a Dirac Hamiltonian with a mass ternf^. 
By calculating the energy spectrum of this Hamiltonian, 
we can relate to the bands corresponding to k l . In 
Fig. [3] the four energy bands are depicted inside and out- 
side the barrier region. Bands belonging to the same 4*4^ 
have the same colour. From this we can see that within 
different energy ranges the transmission at normal inci- 
dence depends on the availability of states corresponding 
to the same I. 

In the two band approximation, and neglecting the 
skew hopping parameters, a similar symmetry transfor- 
mation leads to the Hamiltonian 
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For normal incidence and without bias, the system cor- 
responds to that of two non interacting Schrodinger par- 
ticles described by the Hamiltonian 
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and the wave functions 
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where I = ±1 and m = ~{i/2v 2 F . Since the system is now 
described by a two component spinor one can extend 
the spinorial analogy of the charge carriers in monolayer 
graphene by introducing a pseudospir^. The value of I 
now corresponds to the pseudospin state of the particle. 
As shown above the pseudospin is the consequence of 
the (anti)symmetric nature of the wave functions under 
in-plane mirroring. Although the analogy with the nor- 
mal spin properties of an electron are more pronounced 
within the two band model, the four band model incor- 
porates the same symmetry and therefore the notion of 
pseudospin is applicable in this model as the value of I. 
Note that these symmetry considerations also hold when 
the skew hopping parameters are included because the 
Hamiltonian of Eq. (14) is block diagonal at normal in- 



cidence for a system without bias. 

The different energy ranges in which at normal inci- 
dence and for (5 = new phenomena are expected are 
defined by the value of e = (E — V)/hvp and thus by 
the difference between the kinetic energy and the height 
of the potential in the region under consideration. When 
V — Ti < si < Ti outside the barrier region, but en < 
inside, for positive £/ the state which propagates out- 
side will not be able to propagate inside the barrier. The 
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FIG. 3. (Colour online) Schematic representation of the en- 
ergy spectrum inside and outside the barrier region. The dots 
and arrows indicate the energy regions for which electrons im- 
pinging perpendicularly on the barrier will be transmitted or 
reflected. 



absence of propagating ty + states inside the barrier sup- 
presses the transmission in this energy region. The fact 
that the transmission is suppressed even though there are 
propagating states inside the barrier has been noticed 
before and was called a cloaking of the , F_ states^ and 
is present in both models. 

When Ej > Ti outside the barrier, the states can 
also propagate in regions I and III. Inside the barrier 
however, these states are trapped and their energy spec- 
trum is discrete. This means that when electrons prop- 
agate through the , F_ state, they can propagate inside 
the barrier region only if the energy matches one of the 
discrete energy levels inside. The discretisation condition 
corresponds to Fabry- Perot resonances, namely d = n\/2 
where A = 27r/fc ± . This results in resonances whose en- 
ergies are, in the four band model without skew hopping 
or interlayer bias, given by^ 
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The above reasoning is valid for (near) normal incidence. 
The Hamiltonian of Eq. ( 14 ) however mixes up the two 



states for non-normal incidence. 

To calculate the transmission probabilities at non- 
normal incidence, the transfer matrix method together 
with appropriate boundary conditions was implemented. 
The plane wave solution for the Schrodingcr equation of 
the four band model is a four component spinor and can 
be represented by a product of matrices 



* = VSC, 



(21) 



in which £ corresponds to a 4 x 4 diagonal matrix con- 
sisting of exponential terms and V is for 6 = given 
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in which kr^ is defined in Eq. Q. To find the transmis- 
sion and reflection probabilities, one has to equate the 
wave functions at the borders of the potential barrier. 
This results in two times four equations from which we 
obtain the components of the vector C. Considering the 
boundary conditions of the system, this vector is given 
by 
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where / indicates the wave vector ^ and 5z,±i is the Kro- 
necker delta. Using the transfer matrix method finding 
the coefficients in these vectors corresponds to solving 
the matrix equation 



Cj — Mi^iiMji^ijiCjjj. 



(24) 



Mi^n is the transfer matrix from the region before the 
barrier to within the barrier and Mjj^jjj is the one from 
within to behind the barrier. Using the matrix form of 



the spinor wave functions in Eq. ( 21 1 , the transfer matrix 
is given by 



(25) 



Finally the transmission (T) and reflection (R) probabil- 
ities are obtained as 
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This takes into account the change in velocity of the 
waves when they are scattered into a different propaga- 
tion mode. 

Some of the different probabilities in Eq. ( 26 ) can be 



related to each other via the time reversal symmetry of 
the system. The Hamiltonian under consideration de- 
scribes electrons in the vicinity of one of the two Dirac 
points in reciprocal space which are called the two val- 
leys. The Hamiltonian H' for the other Dirac point is 
given by 



H' = -H 7 



(27) 



This means that electrons scattering from k + to k~ when 
moving from left to right in the first valley arc equivalent 
to electrons scattering from kr to k + but moving in the 
opposite direction and in the other valley. When 6 = 0, 
both valleys are equivalent and time reversal symmetry 
holds near a single Dirac point. Therefore, the trans- 
mission probability of electrons moving in the opposite 
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direction must be the same because of the valley equiva- 
lence. Using a similar argument for the reflection in both 
valleys, one can conclude that 
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Another symmetry operation ensures the symmetry of 
the probabilities with respect to normal incidence. Note 
that the Hamiltonian of the system given in Eq. ([lj is 
not symmetric under a sign flip in k y . An interchange of 
the Ai and B2 atoms together with B\ <H> A2 however 
leaves the system invariant, but corresponds to exchang- 
ing k y — > — k y . Since the system is invariant under this 
transformation, he transmission and reflection should be 
symmetric with respect to normal incidence. The appli- 
cation of an interlayer bias however breaks this exchange 
symmetry and therefore asymmetric results in the trans- 
mission and reflection are expected. This asymmetry was 
noted before by Nilsson el a/pS. The application of an in- 
terlayer bias furthermore lifts the valley degeneracy and 
with it the above discussed symmetry in the scattered 
transmission. Note that when the same quantities are 
calculated for states in the other valley, the asymmetry 
is reversed and therefore the overall symmetry of the sys- 
tem is preserved. 

In the two band model, one obtains a two spinor which 
can be described by a similar matrix product as before 
with 
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where k and k are defines as before. Equating both wave 
functions and their derivatives at the edges of the barrier 
leads to the transmission probability T. 



IV. NUMERICAL RESULTS 

In order to investigate the importance of the skew hop- 
ping parameters, we present in Fig. [4] the transmission 
probability through a potential barrier at normal inci- 
dence with and without including the skew parameters 
as function of the Fermi energy and the width of the bar- 
rier. The results show that even in the low energy range 
where the effect of the skew hopping parameters are ex- 
pected to be the largest, the transmission probabilities 
are very similar. We therefore conclude that for the dis- 
cussion at hand these parameters are not important and 
we will neglect them in the following discussion. 

In Fig. [5]we show the transmission and reflection prob- 
abilities for a pn-j unction as a function of the energy of 
the incident wave E and its transverse wave vector k y . 
The height of the potential is set to Vq — I71 and the 
interlayer potential difference S = 0. The results show 
qualitatively different regions in the (k y , -E)-plane which 
can be explained by identifying which modes are propa- 
gating inside and outside the pn-junction. The borders 
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FIG. 4. (Colour online) Comparison of the transmission prob- 
ability through a barrier of height Vo = O.O571 with (dashed 
curves) and without (solid curves) the skew hopping parame- 
ters. (Left) The energy dependence of the transmission prob- 
ability for junctions of width d = lOmn (blue), d — 25nm 
(red) and d = fOOnm (green). (Right) The width dependence 
of the transmission probability for a Fermi energy of E — | Vo 
(blue), E = §Vb (red) and E = f V . 



between these regions are indicated by dashed curves su- 
perimposed on the density plots. 

For normal incidence, k y — 0, the expected cloaking in 
the and channels occurs for Vo — 71 < E < Vo- 
When the transverse wave vector k y differs from zero, 
it is possible for electrons to scatter into the k~ propa- 
gating mode inside the pn-junction and this results in a 
scattered transmission in the Tj~ channel. For energies 
smaller than Vq — 71, there are propagating /c + -states in 
the junction and this leads to a non zero transmission in 
the Tt channel. This transmission is absent when the po- 
tential Vq < 71 and it is also not present in the two band 
model. For energies larger than the height of the bar- 
rier, E > Vo, the particles behave similar to Schrodinger 
particles. 

The numerical results of the probabilities T~ and R~ 
arc shown in the fourth row of Fig. [5] For E < Vq the 
electrons tunnel using the propagating k~~ states in the 
junction since the states match the number I = — 1. For 
Vo < E < Vo + 71 however there are no available k~ 
states and the transmission is suppressed even though 
the energy is larger than the height of the potential bar- 
rier. This is the equivalent phenomenon of the cloaking 
discussed earlier. This phenomenon is absent in the two 
band approximation. 

The scattered reflection (i?^,i?T^) shown in the sec- 
ond and third row of Fig. [5] is very strong at near normal 
incidence and E « Vq. In this region both states are 
propagating outside the junction, but are evanescent in- 
side it. In this situation scattering is favorable due to 
the symmetry of the scattering probabilities. This ar- 
gument is independent of the mode the incident wave 
is in. For a pn-junction the symmetry argument of the 
previous section only holds for the scattered reflection. 
Because inside the junction the carriers behave as hole 
states while outside they are electrons, the symmetry de- 
scribed in Sec. |III|is broken and therefore the equivalence 
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FIG. 5. (Colour online) The transmission and reflection probabilities through a pn-junction of height Vo = §7i as function of 
energy and transverse wave vector k y . The energy is expressed in units of 71 and the wave vector in units of a -1 , the inverse 
of the in plane interatomic distance. The dashed curves indicate the borders between the regions where different modes are 
propagating or are evanescent inside or outside the junction. 



in scattered transmission is no longer valid. Since under 
reflection the electrons return again in an electron state, 
for the reflection channel the symmetry remains valid, 
which is also seen in the calculations of the reflection 
channels i£t and RZ. 

In Fig. [6]we show the transmission and reflection prob- 
abilities for a potential barrier of width 25nm and the 
same height as the pn-junction. The cloaking in the 
and TZ channel occurs for the same conditions as the 
pn-junction. When the transverse wave vector k y differs 
from zero, it is possible for electrons to scatter into the 
k~ propagating mode inside the barrier and this results 
in the observation of resonances which follow the expres- 



sion given by Eq. (20) for k y large. For energies smaller 
than Vq — 71, there are propagating /c + -states in the bar- 
rier and resonances appear which follow the expression 
given by Eq. (20). These resonances are absent when 



the potential Vb < 71 and are also not present in the two 
band model. The symmetry argument of the previous 
section does hold for a barrier since at both sides of the 
barrier the particles have positive energy. Therefore the 
scattered transmission is the same as shown in the second 
row of Fig. [6j 



In Fig. [7] the transmission and reflection probabili- 
is for a barrier of height Vq = §71 and interlayer bias 
= O.371 are plotted. The bandgap introduced by the 
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FIG. 6. (Colour online) The same as in Fig. [5] but now for the transmission and reflection probabilities through a potential 
barrier of height Vo = §71 and width d = 25nm. 



interlayer bias suppresses the transmission in the energy 
region between Vq ± S. In the bandgap a remarkable 
asymmetric feature with respect to normal incidence in 
the reflection channels shows up. This is a manifestation 
of the breaking of the interlayer sublattice equivalence 
as discussed in Sec. |III| This asymmetry is also present 
in the scattered transmission where it depends on the 
incident mode. This asymmetric feature is only present 
when scattering between different propagation modes is 
possible. It is therefore an other qualitative feature that 
is not present in the two band approximation. 

The interlayer potential difference furthermore couples 
the two propagation modes at normal incidence. There- 
fore the suppression due to cloaking is also adjusted. Now 
both the interlayer bias term S as the transverse momen- 
tum k y causes the two modes to be coupled as shown 
in Eq. (14). Due to the interplay of both parameters, 



cloaking at normal incidence splits into two branches at 
finite k y as shown in the non scattered transmission and 
reflection. 

The conductance is shown in Fig. [8] for both the two 
band approximation and the four band method for bar- 



riers of different width (|8^a) and|8^b)) and a pn-junction 
(j8^c) ) . For energies smaller than the barrier's height, 
E < Vo, the resonances in the transmission show up as 
peaks in the conductance. Using the four band method 
those resonances are however more pronounced and there 
is a difference in the position and number of peaks. Fur- 
thermore for energies smaller than Vo — 71 , the resonance 
peaks for the propagation via the fc + -states appear as 
shoulders of the other peaks. This phenomenon does not 
occur in the two band approximation for which all the 
resonances are similar. When the energy is larger than 
71, additional peaks result from propagation via the k~- 
states inside the barrier while peaks of the two band ap- 
proximation do not differ. This is clarified in the inset 
of Fig. showing the contributions of the different 

transmission channels. While the contribution of the 
channel is low in this region, the TZ and the scattered 
transmission channels have well pronounced resonances. 
When the energy is larger than the height of the barrier, 
E > V , the conductance predicted by both models is of 
the same order of magnitude. However, when the energy 
is larger than Vq + 71, the A:~-state is not cloaked any- 
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more resulting in additional conductance that is absent 
for the two band model. 

Snyman et aZp^ have shown that both models coincide 
in the region for low barrier height, i.e. Vq < j±, which 
is confirmed by our calculations. The characteristic form 
of the conductance resembles that of the result obtained 
here for the two band approximation. Now it is however 
clear that the conductance peak just above the barrier 
height is lower in the four band model. This is due to 
the additional scattered reflection channel Rli/R+ that 
is absent for low barriers. In the latter case only the Tf 
channel could contribute and coincides with the two band 
model transmission, but now this channel is suppressed. 

Fig. |8^b) shows that the resonant peaks in the con- 
ductance depend on the width of the barrier but note 
that they occur at different energies in both models. In 
Fig. [8jc), the conductance of a PN junction is calcu- 
lated. Although the barrier is of infinite width, both 
models predict a finite conductance. For the four band 



model however, the different energy regions show up as 
bumps in the conductance caused by the availability of 
the second mode of propagation. 

Fig. j8^cl) shows the conductance for a biased potential 
barrier. The results are similar to that of the unbiased 
case but are influenced by the suppression of the conduc- 
tance in the energy range of the bandgap at Vq ± S. Note 
that although the scattered transmission is asymmetric 
with respect to normal incidence, this is not visible in the 
conductance calculations since it sums over all values of 



V. CONCLUSION 

We evaluated the transmission and reflection of elec- 
trons through potential barriers and a pn-junction in bi- 
layer graphene. We extended previous calculations per- 
formed within the two band model to the four band 
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model. We compared the results with and without 
taking into account the skew hopping parameters and 
found that the latter can be neglected for energy ranges 
E > 71/IOO « 4meV. Within the four band model, the 
results from the two band approximation are recovered 
for small energies and low potential barriers. We find 
new phenomena such as transmission resonances at nor- 
mal incidence. We showed that the notion of pseudospin 
used to describe electrons in the two band approximation 
corresponds to the wavefunctions being symmetric or an- 
tisymmetric with respect to in-plane mirroring and that 
this leads to the observation of cloaking which occurs 
both for symmetric and antisymmetric states. We have 
also pointed out that because it is a consequence of the 
symmetry of the system, the notion of pseudospin also 
holds when the skew hopping parameters are taken into 
account. For high energies a new mode of propagation is 
available for the electrons, which is not present in the two 
band approximation, and we found that for non normal 
incidence it is possible to scatter between the two modes. 
The resulting conductance of the four band model incor- 
porates these new phenomena and therefore differs signif- 
icantly from the conductance calculated within the two 
band model. This difference manifests itself by the pres- 
ence of many more and well defined resonances and a sub- 
stantially higher conductance for high energies. Finally, 
we showed that the application of an interlayer bias in 
the system significantly changes the transmission charac- 
teristics. The bandgap created by the interlayer poten- 
tial forms a distinct feature in the conductance and the 
transmission and reflection probabilities lack reflection 
symmetry due to normal incidence. The introduction of 
the interlayer symmetry breaking term furthermore cou- 
ples the symmetric and antisymmetric modes. Therefore 
the notion of pseudospin is no longer valid. 
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FIG. 8. (Colour online) Energy dependence of the conduc- 
tance of a single barrier with d — 25nm (a), d = Wnm (b) 
and a pn-j unction with height Vo = §71 calculated using the 
four band method (blue) and the two band approximation 
(dark red). Figure (c), (d) and the inset of (a) show the 
different contributions of the four transmission channels as 
dashed lines, (d) Conductance for a biassed bilayer with in- 
terlayer bias S = O.371 and potential barrier the same width 
and height as that of (a). 
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